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Abstract 

The Villain form of SO(3) lattice gauge theory is studied and compared to Wilson’s 
SU(2) theory. The topological invariants in SO(3) which correspond to twisted 
boundary conditions in SU(2) are discussed and lattice observables are introduced 
for them. An apparent SO(3) phase with negative adjoint Polyakov loop is explained 
in terms of these observables. The electric twist free energy, an order parameter for 
the confinement-deconfinement transition, is measured in both theories to calibrate 
the temperature. The results indicate that lattices with about 700 4 sites or larger 
will be needed to study the SO(3) confined phase. Alternative actions are discussed 
and an analytic path connecting SO(3) and SU(2) lattice gauge theory at weak 
coupling is exhibited. The relevance for confinement of the centre of the gauge 
group is discussed. 


1 forcrand@phys.ethz.ch 

2 j ahn@phys. et hz. ch 




1 Motivation 


Although SO(3) lattice gauge theories (LGT) are expected to have the same con¬ 
tinuum limit as theories with gauge group SU(2), they are interesting for at least 
two reasons. First, the equivalence of the non-perturbative continuum limits has 
not been proven yet; there are in fact some indications that they might be different. 
Second, they provide a tool to study the role of the centre of the gauge group, which 
is often considered to be closely related to confinement: the centre of SO(3), as 
opposed to that of SU(2), is trivial. 

The expectations on the continuum limit are largely based on the naive con¬ 
tinuum limit as obtained by expanding the classical action around its minima. As 
SO(3) = SU(2)/Z 2 is locally indistinguishable from SU(2), in particular in any 
small region around the identity, plaquette actions based on SO(3) and SU(2) link 
variables have the same naive continuum limit. Invoking standard arguments of uni¬ 
versality of the continuum limit, one expects the non-perturbative continuum limit 
to be the same, too. However, this expectation is not based on a rigorous theorem 
and should be checked for each individual case. Studies of the Stefan-Boltzmann 
law in both theories for instance have raised doubts jTj. 

Studies of mixed actions containing both SU(2) and SO(3) parts [the latter being 
considered as the adjoint representation of SU(2)] have not yet been able to demon¬ 
strate this universality. Rather they reveal a line of first order phase transitions 
separating the weak coupling regions of the pure SO(3) and SU(2) theories [2, [3j, so 
the continuum limits are not analytically connected in this plane of couplings. The 
location of the phase transitions is essentially independent of the lattice size, so they 
are considered as bulk. In the pure SO(3) theory, the transition occurs at very weak 
coupling. This puts an upper bound on the lattice spacing in the weak coupling 
phase, and prohibitively large lattices are required to study low temperatures. 

A further peculiarity of the SO(3) theory is the occurance of a new meta-stable 
“phase” (long-lived state in the Monte-Carlo simulation) with a negative expectation 
value of the adjoint Polyakov loop at weak coupling nnsi Apart from the Polyakov 
loop, this “phase” seems to be very similar to the standard phase with positive 
adjoint Polyakov loop. Still its existence has fuelled speculations about differences 
between SO(3) and SU(2) theories. We will show that the mysterious phase can be 
understood in terms of the SU(2) theory. 

The question of whether SO(3) and SU(2) LGT describe the same quantum 
theory is also important for our understanding of the role of the centre of the gauge 
group. A common belief asserts that deconfinement is caused by the breakdown 
of the centre symmetry ,fJJ. However, it has been argued that this symmetry is a 
peculiarity of the lattice theory and not present in other regularisations j7j . In the 
SO(3) lattice theory, the centre symmetry is absent and cannot be broken. If this 
theory has a confinement-deconfinement transition, too, it has to be explained by a 
different mechanism. 

While the centre symmetry is absent, it should be noted that centre vortices do 
exist in an SO(3) theory: their definition actually does not rely on the centre but on 
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a non-trivial first homotopy group of the gauge group modulo its centre. As SO( 3) 
is 577(2) modulo its centre Z 2 , this is the same for both groups, 7 Ti[SU(2)/Z2] = 
7Ti[SO(3)] = Z 2 . So confinement mechanisms based on centre vortices can apply to 
both theories equally well. 

In the absence of the centre symmetry, the question of an order parameter for 
the confinement-deconfinement mechanism arises. The Polyakov loop in the fun¬ 
damental representation is useless because its sign is not determined. The adjoint 
Polyakov loop is not good either, since it has a non-zero expectation value also in 
the confined phase because adjoint charges can be screened by gluons. A useful 
alternative is known for the SU(2) theory. The expectation value of a large spatial 
’t Hooft loop follows a perimeter law in the confined phase and an area law in the 
deconfined phase, the coefficient defining the “dual string tension” m- A partic¬ 
ularly convenient way to measure ’t Hooft loops is provided by twisted boundary 
conditions on a hypercubic lattice (a torus). These introduce chromo-electric and 
-magnetic fluxes through the torus and correspond to ’t Hooft loops of maximal 
size EH- Their expectation value has been measured both at zero EU and at finite 
temperatures mm, and perimeter respectively area law have been demonstrated. 
We will show that the free energy of these fluxes can also be measured in an SO(3) 
LGT. As ’t Hooft loops can be thought of as generating centre vortices, this shows 
that the latter do indeed play a role. Preliminary accounts of the present work have 
been given in mm- 

We begin in Sec. |2] with a review of the topological invariants of SU(iV) and 
SU(AQ/Zat lattice gauge theory on a 4-torus, and of the implementation of twisted 
boundary conditions in SU(iV) LGT. Section |3| recalls a formulation of the SU(2) 
theory as a coupled SO(3)-Z2 theory. This formulation motivates a definition of 
twist in the SO(3) theory which is presented in Sec. 0] In Sec. 0 an interpretation 
of the “phase” with negative adjoint Polyakov loop is presented. Section 0 describes 
our numerical computation of the free energy of electric twist in SO(3). Before the 
conclusions, a discussion of alternative actions and an analytic path between SO(3) 
and SU(2) at weak coupling are proposed. The general relation between Polyakov 
loops and twist is reviewed in the appendix. 

2 Topology in SU(iV)/Zjv and SU(AT) on the torus 

In the continuum, there is no local difference between SU(iV) and SU(iV)/Zjv gauge 
fields. The gauge potential can be expressed either in the fundamental representation 
(with generators t a ) or in the adjoint one (with generators T a ), 

Al = Ap? or = A‘T a . (1) 

The components are always the same. The only difference is in the boundary 
conditions, since the transition functions take values in the gauge group. On a 
torus with extent L v in direction e v : 

A l /A ( x + L ^) = ®% /A (. x ) (^ /A 0) - id ») ^ /At 0d • ( 2 ) 
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The transition functions satisfy the consistency conditions 


( 3 ) 

( 4 ) 


ft*(x + L v e v ){ll(x) = z liV Sll(x +Lpe^) to*(x) 

+ L v e v ) fl^(x) = Q^(x + Lfje/j) itf(x) . 

The twist z Ml/ = exp (^pn^) with integer n Ml/ has to be an element of the centre 
Zjv of SU(N). The group SU(IV)/Zjv has trivial centre, so no twist is allowed for 
With © and ©, both descriptions are equivalent also globally, as any set 
of transition functions in SU(IV)/Zjv can be lifted to SU(iV) and any in SU(IV) 
can be projected to SU(N)/Zjy. The twist z = {z^}, specified by the consistency 
conditions in the case of SU (TV) fields, becomes an ordinary topological invariant 
like the instanton number for SU (N)/Jj N fields. In a slight abuse of language, we 
shall call it twist also there. In addition to the 6 twists z )W € 7ri[SU(IV)/Zjv] = 
which are associated with two-dimensional sections of the torus (//^-planes), there 
is one more topological invariant q £ 7T3[SU(iV)/ZW] = Z associated with the bulk 
m mi- The usual topological charge (instanton number) Q is a combination of 
both, 

h\j 

Q = q- — with k = | £^ P an^n pa . (5) 

It can take fractional values for non-zero n. 

On a lattice, periodic link variables represent all possible boundary conditions 
with trivial consistency conditions [i.e. © or m with z /;i , = 1], Therefore, the 
SU(iV)/Zjv lattice theory (with periodic boundary conditions) automatically con¬ 
tains all “twist” and instanton sectors while the periodic SU (N) theory contains all 
instanton sectors but only trivial twist. A non-trivial (fixed) twist z can be imple¬ 
mented by modifying the periodic boundary conditions to U u (x + L^) = z fW U u (x) 
for fi < v for a single value of x^. Alternatively, one can work with periodic link 
variables by absorbing the factor into one plaquette in each pz/-plane (rSJ. The 
partition function of SU(AQ LGT with Wilson action for twist z then becomes 


yWilson 

Z/ SU(AT) 


(*) 


/ TT^xp 

'SU (N) z 




where l labels links, p plaquettes and 


Qx,llV 



if x M = x v 
otherwise. 


0, 


( 6 ) 


( 7 ) 


In the SU(A')/Zat lattice theory, the twist cannot be imposed in this way by the 
boundary conditions - all twists z are included within periodic boundary conditions. 
A “field-theoretic” definition is not possible either, since z does not have a local rep¬ 
resentation as an integral like the topological charge (Q oc f FF). S However, an 

3 In terms of U(iV) fields they do m, but we will not pursue this any further. 
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SU(iV)/ZW lattice theory can also be described redundantly in terms of SU(iV) vari¬ 
ables. Just imagine mapping each link variable to an SU(iV) element, with arbitrary 
Zjv phase. The resulting link configuration will in general not be periodic, but only 
periodic modulo Zjv- The transition functions will satisfy twisted consistency con¬ 
ditions like ©, but the twist will be a priori arbitrary. One way of assigning a 
unique value, would be some smoothness requirement. One could, for instance, try 
to fix the ambiguity by demanding that the SU(iV) plaquette variables be as close 
to 1 as possible. This is not natural, however, since SU(iV)/Zjv lattice gauge fields 
describe continuum gauge fields with Zjv monopoles. Attempting a “smooth” lift to 
SU(iV) ignores the latter. We therefore choose a different approach, defining twist 
by analogy with the SU(iV) theory in a formulation in SU(iV)/Zjv and Zjy variables. 


3 SU(2) LGT as a coupled SO(3)-Z2 theory 

In order to carry the definition of twist over to an SU(iV)/ZW LGT, we invoke a 
formulation of the SU(iV) theory in terms of SU(iV)/Zjv and Zjy variables. We 
specialise to iV = 2, the case studied numerically in this paper. 

The Wilson partition function for twist z © can be expressed as ©a hhe 2 esi 

=v y, I n^,i) n*(«v.v) • 

Ctp = ±l l ' p ' c p<v 

( 8 ) 

where c labels cubes and 

Oi c 11 OL<p i 
pEdc 

with x p fixed for p ^ n,v in the last equation. The factor A f represents a change 
of normalisation. The choice of Z 2 plaquette variables here is different from that of 
|2HE2|: they are related as a p = a p sgn Tr U p . 

Equation © can be understood as follows. The first set of constraints (on a c ) 
can be interpreted as the Bianchi identity for the field strength a p of a Z 2 gauge 
theory. It forbids monopoles and ensures local integrability: a p can be obtained 
from a ZL link field (potential) 7 /, 


( 9 ) 

( 10 ) 


“r = II 9 ' 1 (locally). ( 11 ) 

ledp 

This is possible whenever the product of a p over all closed surfaces is 1. In R 4 , 
all closed surfaces can be constructed from elementary cubes, so the first set of 
constraints suffices. On a torus, there are 6 independent winding surfaces which 
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cannot be obtained in this way. The second set of constraints fixes the products of 
a p on these, ff z pv — 1, a p can then be integrated globally. Otherwise, one has to 
make a pv trivial first. This is achieved by dividing a p by any representative with 
the same a pil , for instance CD; so the potential is introduced as 

a p=c P n 11 ■ ( i2 ) 

l&dp 

With this representation, the sum over a p in © can be replaced by one over 7 p 
There is no Jacobian because the Z 2 gauge orbits all have the same size. The 
constraints are now redundant, so 7 / only appears in the action. This dependence 
can be removed by a change of variables Ui —> 7 ^Ui which replaces a p by the fixed 
(p . Due to the invariance of the Haar measure, the 7 1 sum can then be done trivially, 
and the Wilson partition function © is recovered. 

4 Twist in SO(3) LGT 

The bridge to SO(3) is now built by noting that, without the constraints, © becomes 
the Villain partition function of SO(3) LGT |24j . 

Z so(3j n = X] [ II dUl exp (I q p TrF U p] = / II <iUl exp ( 

ap=±1 Jsv(2) , p ) J 1 

The Villain action, obtained after summation over a p , 

Sy[u ] = - In[2 cosh(f Tip U p )] (14) 

is invariant under Ui —> — Ui for each link separately, so (HTT1) defines an SO(3) LGT, 
redundantly formulated in terms of SU(2) variables. 

A comparison of Eqs. © and cl exhibits the differences between SU(2) and 
SO(3) LGT. The first set of constraints in © shows that the SO(3) theory contains 
local degrees of freedom not present in SU(2): the Z 2 monopoles of the field a p . 
Loosely speaking, these monopoles are responsible for the bulk phase transition 
separating phases with and without monopole condensation in SO(3) LGT. Since 
a p is only an auxiliary field, it is more appropriate to say that the SO(3) monopoles 

Vc = sgT1 TrF U p ( 15 ) 

pGdc 

condense at small (3. In SU(2), these are strongly suppressed because negative 
plaquettes carry a large action. 

The second set of constraints shows that, as opposed to SU(2), the partition 
function of SO(3) LGT with periodic boundary conditions contains all possible topo¬ 
logical sectors, with arbitrary instanton number and twist z. This is as expected 


-S v [t/]) ■ 

(13) 
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because in SO (3) all sectors are represented by transition functions with standard 
consistency conditions Eq. 0- 

These considerations lead us to a natural definition of the twist z in SO(3) on 
the lattice. In the formulation (JHJ) of SU(2), the twist is given by Eq. (fTUl) in terms 
of the auxiliary variables a p . In SO(3), it is better to define it in terms of U pi 

Vpv y j 'y ^ | sgn Trp U X}P ,v (Spi/pa — 1) ■ (16) 

-L J p-*- J (T 

X p ,X G Xp,Xi/ 

The product measures the Z 2 flux of the SO(3) gauge field C//Z 2 through the fiu- 
plane with fixed x p and x a . The apparent dependence on the sign of the SU(2) links 
Ui drops out because each link appears in two plaquettes in the product. Thus r] lw is 
a proper SO(3) observable. Contrary to (fTUl) . it can be defined in any SO(3) lattice 
gauge theory in formulations with or without auxiliary fields a p . The average over 
all parallel planes with a given orientation has been introduced because in SO(3), 
as opposed to SU(2), the flux can change from plane to plane due to the monopoles 
(USD- This means that r] pi/ can take fractional values — 1 < r) pi/ <1. A Z 2 -valued 
twist can be defined by dividing the range of into two pieces: 

Zpviv) = sgn (rjuu, - c pu ) (17) 


with constant c pu . Note that there is no symmetry r] pv —r) pu , so c pi/ = 0 is not 
singled out. The choice of c pv will be discussed in more detail later. The ambiguity in 
assigning a twist z to a given configuration is similar to that of the usual topological 
charge in SU(2). It is caused by lattice artifacts (the monopoles 7/ c ) which make the 
topological sectors connected. What is important is the behaviour in the continuum 
limit (3 —> oo. Since monopoles are exponentially suppressed in this limit, one can 
hope that z pv {rj) becomes independent of c pi/ there. 

Consider first the naive (classical) continuum limit. In sectors with so-called 
orthogonal twist, n = 0 mod N [cf. (JUJ)], the configuration with lowest action (“clas¬ 
sical ground state”) actually has zero action even though twist is present 121 Ell¬ 
in SU(2), this “twist eater” can be constn 

u„(*) = 

where the 4 matrices satisfy 

T T T t T t 

This generates the plaquette field U p = (* which just cancels the factor ( p in the 
Wilson partition function with twist (JUJ), so the action vanishes. The same config¬ 
uration minimises the Villain action m, because the latter is independent of the 
sign of U p . The observable r] pu (fTUl) . evaluated in this configuration, is identical to 
the SU(2) twist z pi/ of (fTUl) . The same holds true for small fluctuations around m , 


Lcted as 

if x p = 0, 
otherwise, 


= z; u i . 


(18) 


(19) 
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so 77 pi/ has the proper naive continuum limit. For non-orthogonal twist, the lowest- 
action configurations are not known analytically, but it is expected that they still 
satisfy sgn Try U p = a p provided the lattice is fine enough. So 7j pi/ again coincides 
with the SU(2) twist. 

Non-perturbatively, i) pi/ is affected by lattice artifacts (rj c ) - just as the geometric 
definition of the instanton number is m- In our Monte-Carlo simulations, the 
distribution of was always divided into two well-separated peaks close to ±1 
(see Fig. [21 further below), so an unambiguous definition of twist was possible. For 
very large volumes, however, a (centre-blind) restriction of the plaquette angle might 
be necessary, like for geometric definitions of topological charges. The Z^-invariant 
choice 

| Try U p \ > e where £ > 0 (20) 

makes the twist sectors disconnected (with c pv = 0 in (fT71) h 4 For sufficiently large 
e, lattice artifacts are completely suppressed and 77 ^ = ±1. In this case, however, 
the modified Villain theory becomes equivalent to an SU(2) theory with plaquette 
restriction Tiy U p > £, summed over all twist sectors. The value of this theory for a 
comparison of SO(3) and SU(2) lattice gauge theory is thus questionable. We shall 
not apply a constraint like (EDI). For the parameters used in our simulations, the 
ambiguity in m has negligible impact. 


5 Negative Polyakov loops 

In this section, we would like to comment on the apparent weak coupling phase 
with negative adjoint Polyakov loop observed in Refs. mm In these articles, SO(3) 
lattice gauge theory was studied with both the Wilson and the Villain action. For 
the lattice sizes used (N t < 8), only a single, first-order phase transition was found. 
This is the bulk phase transition below which monopoles condense. Although its 
position was found to be essentially independent of N t , it was also interpreted as 
a deconfinement transition, because the expectation value of the adjoint Polyakov 
loop is very small below the transition and has a large modulus above. There, two 
different meta-stable “phases” were observed, depending on the initial configuration 
of the Monte-Carlo simulation. The two phases are distinguished by the sign of 
the adjoint Polyakov loop. The positive value can be understood as a standard 
deconfined phase where | Trp P has an expectation value tending to ±1 in the limit 
(5 —> 00 . The adjoint loop Try P = 4(iTiyP) 2 — 1 therefore tends to +3. No 
explanation was found for the other, “mysterious” phase. In the correlator of 
two Polyakov loops and the response to an external field coupled to the Polyakov 
loop were found to be very similar in both phases, so the authors concluded that 
they were physically equivalent. At the critical coupling, tunnellings between both 
weak coupling phases and the unique strong coupling phase were observed. 

4 The discrete auxiliary variables a p are of no concern here, as one can sum over them and work 
with the action <m. 
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Figure 1: Monte Carlo history of the 3 electric twist variables (top) and the adjoint 
Polyakov loop (bottom) (4 4 lattice, [3 = 4.5). Negative Polyakov loop values are 
always accompanied by one or more negative electric twists. 

Based on the comparison between SO(3) and SU(2) LGT, we propose an expla¬ 
nation of the two weak coupling phases in terms of different twist. Consider the 
semi-classical limit. Without twist, all values of the Polyakov loop are compatible 
with zero classical action. These are the “toron” flat directions. The effective po¬ 
tential localises the fundamental Polyakov loop around two values which tend to ±1 
in the continuum limit [3 —> oo. Either way, this implies a positive adjoint Polyakov 
loop tending to +3 in the continuum limit. With non-vanishing twist, the flat di¬ 
rections disappear. For orthogonal twist, the configuration with smallest classical 
action is the twist eater m- The (untraced) Polyakov loop of this configuration 
is just To. The simplest solution of (HU) with a single electric twist £03 = — 1 is 
To = iui and T 3 = icr 2 , so in this case, the fundamental and adjoint Polyakov loops 
are Tr F r 0 = 0 and Tr A Y 0 = —1. In general, the twist eating matrices satisfy (a) 
TrATo = —1 if Zoi = —1 for some i, or ( 6 ) TrA To = +3 if Zoi = 0 for all i, but 
Zij 7 ^ 0 for some i, j. The arguments leading to this known result are collected in 
the Appendix. This classical consideration suggests that the phase with negative 
Polyakov loop is one with electric twist. This can be verified numerically. Figure 
U shows the Monte-Carlo history of the 3 electric twist variables r] 0i as well as the 
adjoint Polyakov loop on a 4 4 lattice at [3 — 4.5, just above the bulk phase tran¬ 
sition at (3 ~ 4.45. The system tunnels between phases with positive and negative 
Polyakov loop. The negative phases always have one or more electric twists and the 
positive phases none. The twist observable (HI) has now been studied also with the 
SO(3) Wilson action, and meta-stable phases with different £ 0? : have similarly been 
found [27] . 

We would like to remark that semi-classically, the above argument is not con¬ 
clusive. The twist eater has a Polyakov loop Tr F P = 0, i.e. on the maximum of the 
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one-loop effective potential. Therefore, it might not be the dominant semi-classical 
configuration; a configuration with a TL^ interface (Euclidean “domain wall”), inter¬ 
polating between the minima of the effective potential, might have a lower effective 
action. Outside of the interface, the adjoint Polyakov loop is positive. For large 
volumes, its spatial average therefore is positive, too. If the configuration with an 
interface indeed has lower effective action, the average Polyakov loop will be positive 
even in the presence of electric twist, provided the spatial extent of the lattice is 
sufficiently large compared to the interface thickness. 

6 Twist free energies 

We shall now measure the free energy of various twists, both in SO(3) and in SU(2). 
The results are used to calibrate the SO(3) lattice spacing. We choose the SO(3) 
coupling /3 = 4.5 just above the bulk phase transition at f3 ~ 4.45 in order to find 
the largest lattice spacing available in the weak coupling phase of SO(3). 

The twist free energies F(z) are defined in terms of partition functions with and 
without twist, 

F{z) = ~ ln W)' ^ 

For electric twist (zoi 7^ 1)) F(z) is an order parameter for the confinement-decon¬ 
finement transition mm- for large volumes, it tends to 0 in the confined phase, 
while diverging as F ~ aL 2 with the spatial extent L in the deconfined phase, a 
is the dual string tension and depends on the temperature. In SU(2), the partition 
functions Z(z) and Z(l) are defined by imposing twisted and periodic boundary con¬ 
ditions, respectively. In SO(3), all twists are summed over within periodic boundary 
conditions, and twist is an observable like the usual topological charge in SU(2). The 
partition functions are defined using m 

Z(z) = zJUSMvM) (22) 

Vo 

where Ztot = | Z(z). In order to determine the value c ^ of Eq. (1171) which 

divides trivial and non-trivial twist, we use the “density of states” as a function of 
a single T) IWl 

P»v(x) = (6(7)^ - x)) . (23) 

Figure [21 shows P 0i for a 8 3 x 4 lattice at fJ — 4.5. We adjust c ^ to the minimum 
of P lw , which is suppressed by several orders of magnitude with respect to both 
maxima, so the precise choice of c has negligible impact, and twist can indeed be 
defined unambiguously, as claimed above. This has been confirmed numerically by 
comparing the above choice of c with c = 0. The results agree within statistical 
errors. 

The strong suppression of intermediate values of i] 0i turns out to be a practical 
problem: a simple Metropolis algorithm is not capable of changing the electric twist 
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Figure 2: “Density of states” Pa of an electric twist (8 3 x 4 lattice, /3 = 4.5). 

with sufficient probability on lattices larger than 4 4 sites. To achieve ergoclicity 
across twist sectors, we use a variant of the multicanonical algorithm |28j . applied 
to r ] 0 i instead of the action (a “multitwist” algorithm so to speak). The depletion in 
the density of states as a function of r] 0i is removed by adding a bias AS(i] 0 i,Vo 2 , V 03 ) 
to the action, so configurations are generated with probability j?e~ s ~ AS . The Monte- 
Carlo process then amounts to a free diffusion across twist sectors, with dynamics 
accelerated exponentially. The modification in probability is corrected by an inverse 
factor e AS in the observables. 

The bias AS(r] 0 i, ?/o 2 , V 03 ) can be represented as a 3-dimensional table since r/ 0? ; 
only takes L 2 + 1 discrete values (cf. Eq. (fT(l h This table is constructed iteratively, 
using converged values on small lattices to form starting values on larger ones. Since 
the barrier is expected to grow approximately quadratically in the system size, the 
initial bias is scaled as A Sy = (L' / L) 2 ASl. Such a table is displayed in Fig. 01 for 
an 8 3 x 4 lattice at (3 = 4.5. If AS were a sum of 3 functions of only one variable, 
the iterative determination would be greatly simplified. Alas, Fig. 01 shows that this 
is not the case. Note that, although (5 has been chosen as small as possible (the 
bulk transition to the strongly coupled phase occurs at f3 ~ 4.45), the necessary 
enhancement of the saddle points reaches 10 12 , as seen more clearly in Fig. 01 which 
is a 1-dimensional cut of the same table along its diagonal (from r] 0i = (—1, —1, —1) 
to (1,1,1)). There, the resulting flatness of the Monte Carlo sampling probability 
is clearly visible. The measured density of states shows two well-separated peaks, 
corresponding to twist sectors (1,1,1) (analogous to SU(2 ) pbc ) and (—1, —1, —1) 
(analogous to SU( 2) tbc in all 3 temporal planes). The strong suppression of the 
saddle point confirms that ordinary Monte Carlo sampling would remain hopelessly 
“stuck” in one sector, leading to the “phases” observed in earlier studies. However, 
in spite of the great multicanonical acceleration, the Monte Carlo evolution of the 
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Figure 3: The three-dimensional reweighting table used in our multicanonical Monte 
Carlo. Five sections — AS(rj 0 i, rj 02 ) are shown for fixed r/ 0 3 . No twist is at bottom 
left, twist in all 3 planes at top right. The table enhances the probability of sampling 
the saddles between twist sectors. Its entries e -A5 vary by 12 orders of magnitude 
( 8 3 x 4 lattice, (5 = 4.5). 
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Figure 4: Cut of the three-dimensional reweighting table along its diagonal (red +). 
The result of multicanonical Monte Carlo sampling is a nearly flat histogram (green 
x). The density of states shows the dominant twist-0 and the smaller twist-3 sectors 
(blue *). 
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Figure 5: Comparison of the twist free energies F(z) between SO(3) at (3 — 4.5 and 
SU(2 ) at various (3. Electric twist in 1, 2 and 3 planes is considered, on lattices of 
size 4 3 x 4, 6 3 x 4 and 8 3 x 4. One finds (3so{ 3 ) = 4.5 <—> Psu( 2 ) = 4.13(3). 


twist variables r/ 0? ; is still slow, and simulating a 10 3 x 4 lattice remains beyond the 
edge of our computer resources. 

This prevents us from reaching spatial sizes large enough for a reliable deter¬ 
mination of the dual string tension |H|. Nevertheless, we can still compare the free 
energies for various small volumes with those obtained in SU(2) with the method of 
j9J. This also allows us to calibrate the SO(3) lattice spacing. 

The twisted free energies F(z) are continuum quantities, which depend upon the 
spatial size L and the temperature T. If the SO(3) and 50(2) theories represent 
the same continuum physics, then to each (3so{ 3 ) should correspond a value (3su(2), 
which yields the same lattice spacing and thereby the same F(z) for equal lattice 
sizes, modulo small lattice artifacts. We test this statement in Fig. 0 The 9 free 
energies F(z) for 1, 2 or 3 electric twists measured on 4 3 x 4, 6 3 x 4 and 8 3 x 4 
lattices in 50(3) at (3so( 3 ) = 4.5 are compared with similar quantities measured in 
SU(2) at different values of f3su(2)- The straight lines show a linear interpolation in 
/3su( 2 ) of the 577(2) data. One observes a very good match (x 2 /dof ~ 1.35) of all 9 
observables, for Psu( 2 ) ~ 4.13(3). 

Thus, a _1 (/350( 3 ) = 4.5) ~ a^ 1 {(3su{2) = 4.13), which is about 200 GeV! As 
conventional wisdom asserts, our 50(3) lattice is very hne indeed. It cannot be 
made coarser because of the bulk transition at (3 ~ 4.45. Therefore, to reach low 
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temperatures T < T c and probe the confined phase would require a lattice of size 
0(7OO 4 ), far beyond what is currently achievable. We would like to emphasise that 
the scale of 200 GeV is in no way related to continuum physics. It is merely due 
to the bulk phase transition of the SO(3) lattice gauge theory beyond which lattice 
artifacts dominate, and gives a lower bound on the cutoffs one can use. This value 
can be shifted by suppressing (or enhancing) lattice artifacts (27J. 

Note the closeness of the matched SO(3)-SU(2) bare couplings = 4.5 vs. 
4.13). This should come as no surprise. Lattice perturbation theory is identical 
between the SU(2) Wilson action and the SO (3) Villain action: the difference resides 
in the ct-monopoles, which do not appear in the perturbative expansion. Therefore, 
Aiattice is the same in both theories, and one should expect similar values for the non- 
perturbatively matched /Ts. The a-monopoles disorder the SO( 3) theory slightly, 
which is why the SU(2 ) matching f3 is slightly smaller. 

7 Discussion and conclusions 

We have seen that a Monte-Carlo simulation of the confined phase of SO(3) lattice 
gauge theory with the Villain action is not feasible in the near future because a bulk 
phase transition imposes an upper bound on the lattice spacing. This is expected 
to be true of the Wilson action as well. Since the bulk phase transition is caused 
by Z 2 monopoles, the question arises whether the confined phase can be reached 
with a modified SO(3) action suppressing them. A chemical potential can move the 
transition to stronger coupling and larger lattice spacing. There are three kinds of 
monopoles in the Villain theory, a c , rj c [cf. Eqs. © and m) and a c = a c 7j c , and 
one can introduce chemical potentials for all: 5 

^(/3, A a , A^, Ao-) e In dUi exp [-S (/?, Aq,, A^, A ct )] (24) 

<Ap=d=l l 

Aq; ? A^, Act) — ^ ^ Tip Up Aq, ^ ^ cx c Xp ^ ^ rj c X a ^ ^ cy. c Tj c . (25) 

p c c c 

The parameter Aq, (for X v = X a = 0) interpolates between Villain SO(3) and Wilson 
SU(2) [cf. ©]. It has been studied by Halliday and Schwimmer [3DJ who found that 
the bulk transition is indeed moved to smaller (3 and disappears (/3 = 0) around 
X c a ~ 0.95, see Fig. El The deconfinement transition can therefore be studied on 
small lattices for A Q > a c q inn. However, since the limit X a —> oo reduces to the 
SU(2) theory, it is clear that A Q ^ 0 reintroduces the Z 2 part of SU(2) missing 
in SO(3). For A Q ^ 0, ( 1251 ) is not really an SO(3) action, but an SU(2) action in 
disguise. It is not suited for verifying the equivalence of SO(3) and SU(2). This 
is expected to be true also of the theory with X a ^ 0, since this also affects the 
variables a p such that they cannot be summed over trivially. 

5 The effect of A Q and X a (with X v = 0) has been studied in |25j . 
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twist 


Figure 6: Phase diagram of Villain SO(3) with er-monopole suppression, from [31]. 


Note that A Q provides an analytic path between the weak coupling phase of 
Villain SO(3) and SU(2), since no phase transition except the above-mentioned 
monopole condensation transition has been found mm- This is a strong argument 
in favour of the equivalence of the two theories in the continuum limit. It is likely 
that an analytic path can be found in the same way for other SO(3) actions (e.g. 
the Wilson action): any SO(3) theory can be formulated with auxiliary plaquette 
variables a la Villain, and related to an SU(2) theory by a chemical potential for 
a-monopoles. The action of this SU(2) theory is to some extent arbitrary: only 
the centre-symmetric part of the plaquette weight exp[— S P (U P )] is required to be 
equal to that of the SO(3) theory. The anti-symmetric part is only restricted by the 
positivity of the weight. We expect the monopole condensation transition to move 
to small /3 with a suitable choice of antisymmetric part for most theories, so that an 
analytic path will exist. 

The remaining option, \ v ^ 0 (but A Q = A CT = 0) does not reintroduce SU(2) 
degrees of freedom: a p can still be summed over, and one obtains a manifestly 
centre-invariant formulation in terms of SO(3) matrices. X v can also be introduced 
in the Wilson SO(3) action. There, like X a in the Villain theory, it shifts the bulk 
transition to smaller, eventually negative (3 [32j . One expects X p to have the same 
effect in the Villain theory also. However, for A^ —> oo the theory now becomes 
equivalent to an SU(2) positive plaquette model: if r/ c = 1 , the plaquette held r] p 
can be expressed in terms of a link held, r) p = Y\dp^ l ~ & A change of integration 
variables Ui —> 7 1 U 1 makes the plaquette variables trivial, rj p = sgn Tr U p = +1. 
This is just the positive plaquette constraint Tr U p > 0. After summing over a p (still 
A a = 0), we obtain a positive plaquette model with action S = — E P {4 Tr £/p+hi[l-|- 
exp(—/3Tr U p )\ }. This modification of action is not expected to have a significant 

6 In the presence of twist, a representative has to be introduced for each twist, like in Eq. Ha- 
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effect on the deconfinement transition. Therefore, comparing SU(2) with Wilson 
action and SO(3) with action S(/3, 0, A, ; , 0), > 0 is little more than comparing the 

former with a positive plaquette version of itself. For the latter, equivalence with 
Wilson SU(2) seems established jJ3]. Moreover, maintaining ergodicity is even more 
difficult than without monopole suppression because the barriers between different 
twist sectors are even higher. 

We conclude that the equivalence of SO(3) and SU(2) lattice gauge theory cannot 
be directly verified by simulation at the present time. The arguments presented in 
the literature against such an equivalence, however, can be invalidated: an analytic 
path between the theories at weak coupling exists under weak assumptions and the 
phase with negative adjoint Polyakov loop observed in SO(3) turns out to correspond 
to sectors with electric twist in SU(2). The only difference between the two theories 
is that a definite twist sector is selected via the choice of boundary conditions in 
5f/(2), whereas all sectors are automatically summed over within periodic boundary 
conditions in SO (3). Also the absence of the centre symmetry in SO(3) is not a 
problem: it is not needed to characterise the deconfinement transition. The (spatial) 
centre vortex - or electric twist - free energy can be used as an order parameter in 
both theories. At the very high temperatures accessible to simulations to date, the 
centre vortex free energies of SO(3) and SU(2) are numerically compatible, provided 
the scales are adjusted between both theories. 

Finally, one can try to generalise the lesson learnt here about the structure of 
the gauge group G necessary for confinement. 50(3) shows that a non-trivial centre 
is not required (cf. 0). On the other hand, the existence of the dual string tension 
characterising the deconfoied phase arises from that of’t Hooft loops or twist sectors. 
Those in turn follow from the non-trivial first homotopy group vri[G/ Centre(G)], 
which is the same for SU( 2) and 50(3). A non-trivial tti[G/ Centre(G)], or more 
precisely tti[G/Zg], where Z G is the discrete part of the centre of G, is necessary 
for a dual string tension to be defined. The non-trivial elements of this group are 
the centre vortex (or twist) excitations. There are only 3 (simple) non-Abclian Lie 
groups for which 7Ti [G/Zq\ is trivial: G 2 , F 4 and E$. The absence of a dual string 
tension in the corresponding gauge theories prompts speculations about the nature 
of the confinement/deconfinement phase transition. The study of G 2 proposed in 
m will be interesting in this connection. 
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Appendix: Twist and Polyakov loops 

Here, we shall review the relation between twist and Polyakov loops for general twist 
eaters, working out the arguments sketched in .To) . 

Consider a twist eater m given in terms of 4 constant matrices G SU(2) 
which satisfy the relations (m, 

= e^T^ (26) 

with = 0,1. These relations define a group, and the matrices can be con¬ 
sidered a (non-unique) 2-dimensional representation of this group. The simple ar¬ 
gument below depends on the fact that, for SU(2), all such representations are 
irreducible unless n^ u = 0 for all p, u. 

A generalised Polyakov loop in the fundamental representation which winds 
around several directions, successively s M times around the //-direction, can be de¬ 
fined as 


Pf{s) = \ Tr F r(s) (27) 

r(s) = r^r'T'T* 0 . (28) 


It is not necessary to consider other orders of the factors P /y ,, since, by Eq. m, t ^ ie 
order only affects the sign of Py (s). Furthermore, irreducibility implies P^ = ±1 [25], 
so |Pf(s)| depends on only modulo 2 and it is sufficient to consider s M G {0,1}. 

We show that, provided n ^ 0, |Pp(s)|, and thus Pa{s), only depend on n /ti , and 
not on the particular representation P^. More precisely, 


M*)\ 


0 if Y.V n ^ s v 7~ 0 ( m °d 2) 

1 Yju n ^ s i' = 0 (mod 2) . 


(29) 


To see this, note that Eq. ( 1351 ) implies 

r(s')r(s) = e i7rs '- n T(s)r(s') ( 30 ) 

where s' ■ ns = y) s' ll n fll/ s u . Now, ns ^ 0 (mod 2) implies that s' ■ ns — 1 (mod 2) 
for some s' and thus 

mom) = -m>mo. (3i) 

There is a theorem that says that the relation m for two invertible matrices T(s) 
and r(s') implies TrpP(s) = 0 (Lemma 1 of [35] . Sec. 4.3). This proves the first 
part of m- If ns = 0 (mod 2), on the other hand, T(s) commutes with all T(s / ), 
in particular the generators If n ^ 0, irreducibility implies that T(s) = ±1, 
proving the second part. 

It is straightforward to convince oneself that the values of |Pp(s)| for = 0,1 
determine the twist uniquely through Eq. (135|) . see Table [T] 

For vanishing twist, all values of the Polyakov loop have minimal classical action, 
since (ESI) is satisfied by arbitrary diagonal matrices These are the so-called 
“toron” modes 55l) . 
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\m\ 

\k\ 2 

Po 

\p\ 2 

P-M 

M 2 

M ■ I< 

\K\ 2 

P-K 

0 

1 

0 

2 

0 

1 

0 

0 

0 

0 

2 

0 

1 

1 

1 

0 

0 

0 

0 

3 

0 

0 

0 

3 

0 

0 

0 

1 

0 

1 

1 

0 

0 

0 

1 

1 

1 

1 

0 

1 

0 

0 

0 

1 

0 

1 

2 

0 

1 

0 

0 

0 

0 

0 

2 

0 

1 

0 

0 

1 

0 

0 

0 

2 

1 

0 

0 

0 

1 

0 

2 

0 

2 

2 

0 

0 

0 

1 

1 

1 

0 

2 

3 

0 

0 

0 

1 

0 

0 

0 

3 

0 

1 

0 

0 

0 

0 

0 

0 

3 

2 

0 

0 

0 

0 

0 

1 

0 


Tabic 1: Multiple Polyakov loops for twist eaters with twist n tJ = eijk'mk , Uqi = 
ki. Polyakov loops winding around up to 2 cycles of the torus are considered: 
P,,. = |Pf(gJ|, ^ = |P F (e 0 + Si) | and M k = |P F (ei + ej)\ with {i,j,k} = {1,2,3}. 
This information is sufficient to determine the twist up to a permutation of space 
coordinates. The individual components can be obtained from the components of 
P, K , M and P 0 , if desired. Some twists are missing because twist eaters only exist 
for orthogonal twist. 
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